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We analyze the desynchronization bifurcation in the coupled Rossler oscillators. After the bifur- 
cation the coupled oscillators move away from each other with a square root dependence on the 
parameter. We define system transverse Lyapunov exponents and in the desynchronized state one is 
positive while the other is negative implying that one oscillator is trying to fly away while the other 
is holding it. We give a simple model of coupled integrable systems that shows a similar phenomena 
and can be treated as the normal form for the desynchronization bifurcation. We conclude that the 
desynchronization is a pitchfork bifurcation of the transverse manifold. 



In this paper we study desynchronization bifurcation for two chaotic Rossler oscillators which are 
mutually coupled via a;-component. They show an interesting behaviour with increasing coupling strength 
(e). When coupling strength (e) is increased beyond some first critical coupling (Eci), the osciUators 
synchronize. As we keep on increasing coupling strength, these oscillators will remain synchronized for 
some time {sci < s < £02)- When coupling strength exceeds a second critical value (£02) the oscillators 
desynchronize. At this point the largest transverse Lyapunov exponent (TLE) become positive. To 
understand this phenomenon in more details, we define systems' transverse Lyapunov exponents (STLE) 
which are specific to each system. In the synchronized state the STLE and TLE have similar value and all 
are negative. But in the desynchronized state one of the STLEs is positive and another is negative which 
I implies that the perturbation grows about one system while it dies out about the other system, i.e. one 

system is trying to fly away while the other is holding it. We present a simple integrable general model 
with quadratic nonlinearity which shows similar phenomena and the nature of this desynchronization can 
be explored in more details with the help of this model. This model can be regarded as a normal form for 
the desynchronization bifurcation in coupled Rossler systems. We also study the cubic nonlinearity and 
find that in this case both SLTEs are negative. 



I. INTRODUCTION 



CO I Chaotic systems when coupled in some fashion or driven by same external signal, synchronize as the coupling 
. strength increases, and are said to be identically or completely synchronized when the variables of the systems become 
O" equal [iH3. We can observe other types of synchronization, such as phase synchronization lag synchronization 0, 
^-H , generalized synchronization etc. Here, we restrict ourselves to identical synchronization, though some of the 

• results may have more general validity. 
^ \ The condition for identical synchronization, or for simplicity synchronization, can be obtained by linear stability 
T '~j ■ analysis. The phase space of the coupled system can be split into two manifolds, the synchronization manifold and 
KnJ ' the transverse manifold. The synchronization takes place when the transverse Lyapunov exponents (TLEs) become 
^. negative [nl, p. 

" " ' An interesting situation arises when two Rossler oscillators are coupled with each other. There are two critical 
coupling constants, Sci, and £c2- For e < Sci, the oscillators are desynchronized. They are synchronized for Ed < £ < 
£c2 and are again desynchronized for e > £c2- In the range £ci < s < £c2, all the TLEs are negative. While outside 
this range one of the TLEs is positive [ll|. Similar considerations apply when one considers a system of coupled 
Rossler systems on a network. When one couples several identical chaotic systems in an array, the desynchronization 
bifurcation at £c2 can be identified as a short wavelength bifurcation, where the shortest spatial wavelength mode 
becomes unstable Ell • 

The purpose of the present paper is to study and understand the desynchronization bifurcation. We find that for 
£ > £c2, the attractors of the two coupled systems split and start drifting away from each other and the rate of 
drift is proportional to \/£ — £c2- We introduce system transverse Lyapunov exponents (STLEs) and we find that the 
largest STLE for one system becomes positive while that for the other system becomes negative. It implies that when 
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desynchronization bifurcation takes place the perturbation transverse to the synchronization manifold grows about 
one attractor while it dies out about the other attractor. Next, we construct a simple model of coupled integrable 
systems which obeys similar properties. We are able to analytically derive the properties of the desynchronization 
bifurcation in our model. The model can be considered to be the normal form for the desynchronization bifurcation. 
In the model we study both quadratic and cubic nonlinearites and we find that STLEs are useful to distinguish 
between the two nonlinearities. The quadratic nonlinerity gives the results corresponding to the desynchronization 
bifurcation in Rossler system. Form this study we identify the desynchronization bifurcation in coupled Rossler 
systems as a pitchfork bifurcation of the transverse manifold. The present paper is divided in the following sections. 
In section III Al we study linear stability analysis of two n dimensional systems. In section IIIBI we define systems' 
transverse Lyapunov exponents and develop an algorithm to calculate STLE. In section FlI CI we presented numerical 
results on Rossler oscillators. We proposed simple integrable model in section IIIIl which shows similar behaviour and 
show how this model can be treated as a normal form of this desynchronization bifurcation. 



II. DESYNCHRONIZATION BIFURCATION 



We first consider the linear stability analysis of the synchronized state of two coupled dynamical systems. Next, we 
introduce system transverse Lyapunov exponents. Then, these are used to study the desynchronization bifurcation in 
the coupled Rossler systems. 



A. Linear stability analysis of synchronized state of two coupled dynamical systems 

Consider an n-dimensional autonomous dynamical system, 

i^fix), (1) 
and couple this system with an identical dynamical system y, 

y = f{y)+e2T{x-y) (2) 

where, ei and £2 are scalar coupling parameters. F is known as the diffusive coupling matrix. In general, F = 
diag(7o, 71, 7n_i), and defines the components of x and y which are coupled. The synchronization manifold is 
defined hy x — y = where s satisfies Eq. ([!]). Let, S^x and be the deviation of x and y from the synchronized 
solution s. We have 

L = V/(s)e,+£iF(e^-^,) 

iy = V/(s)ey + £2F(^.-^,) (3) 
These two equations can be also be written as [l3| , 

^ = v/(s)e + rcG^ (4) 

where, ^ = {S,x,£,y) and G is the coupling matrix. In this case 

Let, Pk be an eigenvector of with eigenvalue fik] G^Pk — fikPk- Operating Eq. ^ on fj,k and defining (k — £,Pk 
and we can write an equation for (^k as [l3| , 

Cfe = [V/(s) + £//feF]Cfe. (5) 
Here, the matrix G has two eigenvalues fio — and fii = — (si + £2)- Thus, Eq. ([5]) gives the two equations, 

Co = V/(s)Co (6) 

Ci = [v/(s)-Aiir]Ci. (7) 

Here Eqs. ([6]) and ([7]) define motion of small perturbations on the synchronization and transverse manifolds respectively 
and these can be used to obtain the Lyapunov exponents for the two manifolds. The synchronized state will be stable 
when all the transverse perturbations die with time, i.e. when all the transverse Lyapunov exponents are negative. 
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B. System's Transverse Lyapunov Exponents 

We now introduce transverse Lyapunov exponents which are specific to the invidual systems x and y. 
The dynamics of the difference vector z = x — y, is 

i = fix) - fiy) - (ei + £2)rz (8) 

In Eq. ([5]) we can expand f{y) in Taylor's series about the co-ordinate of the first system x or f{x) about the second 
system y. This gives us the following two equations. 

i - (z-V)/(x)-(£i+£2)rz (9) 
i = (z-V)/(y)-(ei+£2)rz (10) 

where we neglect the higher order terms. In the synchronized state, Eqs. ^ and PU)) are identical and give the 
transverse Lyapunov exponents. In the desynchronized state, Eqs. ^ and (|10p in general give different exponents 
and we refer to them as system transverse Lyapunov exponents (STLEs) since they are specific to each system and 
denote the largest of them as A^; and Xy respectively. For the synchronized state A^; = Xy and they are negative. For 
the desynchronized state A^; may not be equal to Xy and tell us about how the difference vector z behaves in the 
neighborhood of the two systems. Note that for the synchronized state these STLEs belong to the actual spectrum 
of Lyapunov exponents of the coupled system, but not for the desynchronized state. 



C. Two coupled Rossler systems 

We now take the specific example of two coupled Rossler oscillators [l3| . Denoting the variables of the two systems 
by X and y the coupled equations are 

xi = -X2 - X3 + e{yi - xi) 

X2 = Xi+ arX2 (11) 
X3 = br + X^lxi — Cr) 

and a similar set of equations for the other system y. Here, we have coupled only the first component, i.e. F = 
diag(l, 0, 0) and we take symmetric coupling, e = ei = £2- 

Fig. [1] shows the variation of the largest transverse Lyapunov exponent {Xmax) with coupling strength for two 
mutually coupled identical Rossler oscillators. As discussed in the introduction there are two critical coupling constants 
Eel and £c2- The synchronized state is stable when eci < £ < £c2- At £ = £c2 the system undergo a desynchronization 
bifurcation. As we see in Fig[l] Xmax is positive when e > £c2, which implies the synchronous state is unstable. To 
understand this phenomena in details we calculate the systems' transverse Lyapunov exponents (Xx and Xy) introduced 
in the previous subsection. 

In Fig. [2J the two largest system transverse Lyapunov exponents, A^, and Ay arc plotted as a function of the 
coupling strength e. As noted before, there are two critical coupling constants, £ci and £c2- At both the critical 
points Xx = Xy = 0. For < £ < £ci, the coupled oscillators are desynchronized. The attractors of the two systems 
overlap and are similar in nature. In this region, Aa; ~ Aj^ and both are mostly positive. For, Sci < £ < £c2, the 
two Rossler oscillators are synchronized. Here, Aa, = Xy and both are negative. For e > £c2, the oscillators become 
desynchronized. Here, the largest STLEs show an interesting behavior. One of STLE becomes positive but the other 
becomes negative. Note that for very large values of e the coupled system becomes unstable. 

To understand the result that one STLE is positive and the other is negative, let us first look at the phase space 
plots of the attractors of the two coupled oscillators in Fig. [3^. The two attractors are identical and overlap at £ = £c2 . 
As £ increases the two attractors split and start moving away from each other as shown in Fig. [3^. Figure shows 
the distance D, between the centers of the two attractors as a function of e. For e > £c2, the distance D shows a 
power law behavior, 

^-7(£-£c2)^ (12) 

The fit is shown in Fig. |3}d and the exponent is u ^ 0.474 ± 0.054 ^ 0.5 and the other parameters are 7 = 67.38 ± 
15.75, £c2 = 3.002 ± 0.000004. The power law behavior is a characteristic feature of a second order phase transition. 

Let us now come back to the result of Fig. [2J that for e > £c2 one of the STLE is positive and the other is negative. 
These STLEs tell us about the behavior of the distance between the attractors as viewed from each of them in the 
linear approximation. Thus, we can say that in the linear approximation one of the attractors is trying to fly away 
while the other one is trying to hold them together. The stability of the coupled system implies that the negative 
STLE wins the battle. It appears that as e increases, the hold of the negative STLE decreases and hence the two 
attractors start drifting away from each other and for large values of £ the system becomes unstable. 
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FIG. 1: The largest transverse Lyapunov exponents, Xmax, of the two coupled chaotic identical Rossler oscillators are plotted 
with the coupling parameter e. There are two critical couplings £ci(~ 0.1) and £c2(~ 3.0). In the range Ed < e < £c2 the 
synchronized state is stable. The desynchronization bifurcation takes place when e > ec2- The Xmax is calculated from Eq. [7] 
Rossler parameters are = 0.15, for = 0.2 and Cr = 10.0. Note that for very large couplings the coupled system become 
unstable. 



III. MODEL SYSTEM 



Since Rossler oscillators are chaotic it is not easy to decipher the behavior of the desynchronization bifurcation. 
Hence, we now propose a simple model of coupled integrable systems showing a similar desynchronization bifurcation. 
It is easy to see that with one dimensional systems we do not get the desynchronization bifurcation. Hence the 
minimum dimension is two. The proposed model is 

xi = axi + hx2 + s{yi — xi) 

X2 = cxi+ dx2+ g{xi,X2) 

yi = ayi + by2 + e{xi - yi) 

m = cyi+dy2+ g{yi,y2), (13) 

Here, a,b,c,d,a, l3 are parameters of the systems and g is a nonlinear function of its arguments. As in the case of 
Rossler systems we couple the xi component. The model system is chosen so that the synchronizd state corresponds 
to the fixed point x* = y* = (0,0) for small values of e and we observe a desynchronization transition as e increases 
(isj . For this to happen the parameters of the system must obey the conditions; a + d < 0, d > 0, {ad — be) > 0. 
Under these conditions, the fixed point (0, 0) becomes unstable at the critical coupling constant Ec ~ £c2 — \{o-^ 
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FIG. 2: The largest system transverse Lyapunov exponents, \x,y, of the two coupled chaotic identical Rossler oscillators are 
plotted with the coupling parameter e. The desynchronization bifurcation is observed for large coupling (here at e ~ 3). The 
inset shows a blowup of \x and \y just after the desynchronization takes place. Rossler parameters are = 0.15, br = 0.2 and 
Cr = 10.0. 



A. Quadratic nonlinearity 

We first consider quadratic nonlinearity [l6| . 

g{ui,U2) = a(u\ + f3uiU2 + u\) (14) 

With quadratic nonlinearity, the model system has three fixed points. One is (0, 0, 0, 0) which is also a fixed point 
(0, 0) of the uncoupled systems. The other two fixed points are given by 

[Xi,x2,yi,y2) = (— ±VB, —x^~-y^,A-x^, 

d — E £ 

-here A = - ^(..X^et+a^l-^e)) ^ ^ie) = ^ - ^ - M^A, W ^ + + - Aae - /36(a - 2e). 

For e < Ec, the fixed point (0,0) is stable and it becomes unstable at the critical coupling constant Ec — £c2- For 
E > Ec the coupled system has two stable fixed points given by Eqs. (ITSI) . 

The STLE can be obtained by writing equations for the difference vector z = a; — y as in Eqs. (IH]) and (ITUl) . The 
largest STLEs are given by 

X^^y « ±CV^^c + 0{e - Ec): (16) 
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FIG. 3: a. The projection of the attractors of the two Rossler oscillators (red and blue online) on x-y plane for 
e — 3.000,3.025,3.050,3.100,3.150,3.200. Note that Eqs. pil) obey the x ^ y exchange symmetry and hence the attrac- 
tor obtained by the exchange x <^ y, is also a solution, b. The distance D, between the centers of the attractors of the two 
Rossler oscillators is plotted as a function of the coupling strength e. The continuous curve (blue online) is a power law fit 
(Eq. mi) with the exponent v = 0.5. Note that the diastance D is also proportional to the two solutions obtained by the 
X <^ y exchange symmetry. 



Where C = („^/,^^)^„,.,^-..)^^„(^,_,^^,^^)^ Vf „ F = J^^ > , G ^ d{a + b - (Sab) + (ad - bc)if3b - 2a) 

and H = d^(a^ + 6^ — fiab) — [ad — bc){ad + be — 2/3bd). Figure shows the largest STLE X^^y as a function of the 
coupling constant s. For e < Ec, X^^y are negative and equal. At e = Sc, they are zero and for s > Ec, one of the STLE 
is positive while the other is negative. This behavior of Xx.y is similar to that of the desynchronization transition in 
the coupled Rossler system seen in Fig. [2] 

The distance between the attractors of the two systems, i.e. between {xl,X2) and (2/1,2/2); is given by 

vMEZFE^ + o(.-..) (17) 
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FIG. 4: a. Two largest system transverse Lyapunov exponents, \x,y of the model, Eq. (|13|l with quadratic nonlinearity, Eq. (|14l 
are shown as a function of the coupling parameter e. Other parameters are a = —1.00,6 = — 2.00,c = 1.00, d = 0.50, a = 
— 1.00,/3 = 2.00. b. The distance D between the fixed points of the two systems (Eq. (|15p ) as a function of the coupling 
constant e. 



Figure HId plots the distance D as a function of the coupling constant e. Thus, for e > Ec, D oc ^Je — Sc- 

At the desynchronization bifurcation in the model system the fixed point (0, 0) becomes unstable and two new 
stable fixed points emerge. The distance between the stable fixedpoints grows proportional to — e^- These are the 
characteristic features of the supercritical pitchfork bifurcation [13, [13 • This bifurcation takes place in the transverse 
manifold. This can be seen by noting that the three fixed points of the the model system, can also be obtained from 
the equation satisfied by the transverse component z\ as 

zl{B{E)-{zXf)=^ (18) 

This is a cubic equation and B[e) oc (e — £c) since B{ec) — 0. This equation is exactly the normal form of a pitchfork 
bifurcation .17., 18J . Similar equation can be written for . 

The proposed model with quadratic nonlinearity will show supercritical pitchfork bifurcation when /3 >= 
(ad-bc){a '^~^°'"^) ^ Otherwise it will undergo subcritical pitchfork bifurcation. 



B. Cubic nonlinearity 

We now consider cubic nonlinearity 

g{ui,U2) = a{ul + I3iulu2 + huiu^ + 2:2) (19) 

In Fig. [5] the largest transverse Lyapunov exponent {\max) of this system is plotted with the coupling strength. As e 
crosses the critical coupling strength (ec) the largest transverse Lyapunov exponent become positive and synchronized 
state become unstable. In Fig. [6^ we plot the two largest systems' transverse Lyapunov exponents (A^: and Aj,) of 
the model system given by Eq. [T^lwith cubic nonlinearity as a function of the coupling strength e. Here we can find 
that the exponents have same value for all coupling strengths and everywhere they are negative, except at the critical 
coupling strength, Sc where both of them are zero. 

In the desynchronized state one can calculate the stable solutions analyticallly for the cubic nonlinearity (Eqs. 
and (dH)). The fixed points are given by. 



^&y2d(£-£,) 



_ (a-2e) 

y*2 = -x*2, (20) 
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FIG. 5: The largest transverse Lyapunov exponent, Xmax is plotted with the coupling strength e for the model, Eq. p3p with 
cubic nonlinearity (Eg. llQp . The critical coupling strength is Ec = 1-5. When e > Ec the Xmax is positive and synchronous state 
become unstable. The system parameters are a' = -1.00, b' = -2.00, c' = 1.00, d' = 0.50, a = -1.00, /3i = 3.00, ^2 = 3.00. 



where ec = ^(^ — ^) and F — a{{a — 2e)^ — I32b{a — 2e)^ + ^i6^(a — 2e) — b^}. In the synchronized state the systems 
synchronize in the (0, 0) solution. When the coupling strength e crosses the critical value the systems undergo 
desynchronizaion bifurcation as depicted in Fig. [SJ but all STLEs are negative (Fig. So, the individual systems 
are stable. The distance between the two fixed points is proportional to y/e — Sc and is shown in Fig. ^jp. 
We can calculate the STLEs for cubic nonlinearity by considering the transverse component z = x — y, 



2Fd{F - (a - 2e)G - bH) 
"^'^ ~ {{a + d-2e)F+{e-e,)Gy^"^'' 

higher order (21) 



where, G = 3(a - 2ef + Pib'^ - 2^36(0 - 2e) and 

H = 3b^ + /32(a — 2e)^ — 2/3i6(a — 2e). The STLEs for cubic nonlinearity have linear dependence on the parameter 
after the desynchronization bifurcation takes place and both are negative. 



C. Comparison with coupled Rossler systems 

We now compare the results for the model system with that of two coupled Rossler systems. Comparing Figures 
and|B]3, we see that for both the model and the coupled Rossler systems, for e > £c, D oc — Cc-We note that 
D may be taken as the distance between the attractors of the two systems or the distance between the two solutions 
obtained by the x ^ y exchange symmetry. For the coupled Rossler systems these solutions are chaotic while for the 
model system they are fixed points. The nature of these solutions depends on the synchronization manifold. However, 
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FIG. 6: a. Two largest system transverse Lyapunov exponents, X^^y of the model, Eq. (|13|) with cubic nonlinearity, Eq. (|19|l . 
are shown as a function of the coupling parameter e. Other parameters are a' = —1.00,6' = — 2.00, c' = 1.00, d' = 0.50, a = 
— 1.00, /3i = 3.00, ^2 = 3.00. b. The distance D between the fixed points of the two systems for cubic nonlinearity as a function 
of the coupling constant e. 

the desynchronization bifurcation takes place in the transverse manifold were the both the coupled Rossler systems 
and our model show a very similar behavior. 

For the coupled Rossler systems we can carry out an approximate analysis. We write equations for the difference 
and sum of the variables of the two systems, z = u^^'^ — w*^^^ and s — u'^^^ + u'^^\ and then treat z and s as constants 
near the desynchronization bifurcation. This gives a cubic equation for the transverse components as Z2(/B — zf) where 
B depends on the parameters. The condition B — gives £c2 ^ 3.33 . . . which is somewhat larger than the observed 
value of 3.002 of the desynchronization bifurcation. 

Thus both the transitions in the coupled Rossler systems and our model can be identified as supercritical pitch-fork 
bifurcations of the transverse manifold. 

The nature of the nonlinearity can be identified using STLEs defined by us. Comparing the behavior of STLEs for 
£ > Ec in Figs. [2113^ and [5^, we see that the behavior of SLTEs for the coupled Rossler systems matches with that of 
our model with quadratic nonlinearity, but not with the cubic nonlinearity. 

We find the the form used in Eq. (|T3l) with quadratic (Eq. (|T4l) ) or cubic (Eq. ^T9\\ ) nonlinearity, is the simplest 
form we could get for the desynchronization bifurcation of the transverse manifold. The model also gives the standard 
normal form (Eq. (jl8p ). of the pitchfork bifurcation for the transverse component. Hence, the model may be treated as 
the normal form for the desynchronization bifurcation [l^ . We note that the coupled Rossler systems and the model 
with quadratic nonlinearity have similar properties. Hence, we conjecture that our model of Eq. (jl3p with quadratic 
nonlinearity (Eq. (jl4p ) is the normal form for the desynchronization bifurcation of the coupled Rossler systems. 

IV. DISCUSSION 

From the discussion above, we conclude that the desynchronization bifurcation of the coupled model system, Eq. (jl3p 
as well as the coupled Rossler systems, Eq. (|lip. are supercritical pitchfork bifurcations of the transverse manifold. 
The synchronization manifold decides the nature of the attractor which is chaotic for the coupled Rossler systems 
while it is periodic (fixed points) for our model system. 

We have presented the analysis for symmetric coupling with e = ei £2- If instead we take asymmetric coupling 
£1 7^ £2) the nature of the desynchronisation bifurcation does not change. This is because this bifurcation takes place 
in the transverse manifold defined by the difference vector z and in the equation for z, (Eq. ([8])), we only have the 
sum £1 + £2. We also note that for £1 7^ £2, the x y exchange symmetry exists in the trasverse component though 
not in the longitudinal component. 

We find the the form used in Eq. (jl3p is the simplest form we could get for the desynchronization bifurcation and 
also, we get the standard normal form (Eq. (fT8|) '). of the pitchfork bifurcation for the transverse component. Hence, 
the model may be treated as the normal form for the desynchronization bifurcation. We can further simplify the 
model by choosing a = — l,c=l,a = ±l. We note that the coupled Rossler systems and the model have similar 
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properties. Hence, we conjecture that Eq. (1131) with the quadratic nonhnearirty fEa. (|14p ) is the normal form for the 
desynchronization bifurcation of the coupled Rossler systems. 

Let us now consider coupled Rossler systems on a network. Consider n coupled Rossler oscillators. Denoting the 
variables by u^^\ j — 1,2, . . . , N , the equations can be written as 

iiU) = /(u(j)) + £^ J,fer(u('=) (22) 

where J is the coupling matrix. The analysis of Pecora and Carrol [2l[ shows that the equations for the transverse 
manifold can be cast into a general form of a master equation and is the same as that for the two coupled systems. 
Thus, the present analysis should be applicable for the desynchronization transition for coupled systems on a network. 
How do the attractors of the different systems split for e > ec2? Consider three mutually coupled Rossler systems. 
We observe an interesting phenomena of symmetry breaking. In this case at the desynchronization bifurcation we 
still get splitting of the attractors into two as in Fig. |3K, with two oscillators on one side and the remaining oscillator 
on the other side. The two oscillators on the same side remain synchronized We find that the distance between 
the center of these oscillators varies with the coupling in the same fashion as in Eq. ([T^ . When four oscillators 
are coupled in a reactangle then this desynchronization bifurcation takes place between two pairs of oscillators. The 
oscillators in the same pair remain synchronized. 

V. CONCLUSION: 

To conclude, we have analysed the desynchronization bifurcation in the coupled Rossler systems. We give a simple 
model of coupled integrable systems which shows a similar phenomena. The model may be treated as the normal form 
for the desynchronization bifurcation. After the desynchronization bifurcation the attractors of the coupled systems 
split into two and start moving away from each other. We define system transverse Lyapunov exponents corresponding 
to the difference vector of the variables of the systems. For e > Sc and quadratic nonlinearity, the STLE for one 
system becomes positive while that for the other system becomes negative. While for e > Ec and cubic nonlinearity, 
the STLEs of both systems are negative. From the analysis of the distance between the two attractors which is 
proportional to — £c, the behavior of SLTEs and the cubic form for the transverse components, we conclude that 
the desynchronization bifurcation in the coupled Rossler systems is a pitchfork bifurcation of the transverse manifold 
and has the normal form of our model with quadratic nonlinearity. 
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